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Abstract For a natural exponential family (NEF), one can associate in a natural way 
two standard families of conjugate priors, one on the natural parameter and the other 
on the mean parameter. These families of conjugate priors have been used to establish 
some remarkable properties and characterization results of the quadratic NEF's. In the 
present paper, we show that for a NEF, we can associate a class of NEF's, and for each 
one of these NEF's, we define a family of conjugate priors on the natural parameter and 
a family of conjugate priors on the mean parameter which are different of the standard 
ones. These families are then used to extend to the Letac-Mora class of real cubic natural 
exponential families the properties and characterization results related to the Bayesian 
. theory established for the quadratic natural exponential families. 

CN | Keywords: natural exponential family, variance function, prior distribution, posterior 

expectation. 
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1 Introduction and preliminaries 



To make clear the motivations of the present paper, we first recall some facts concerning 
the natural exponential families and their variance functions. Our notations are the ones 
used by Letac in [9]. Let \i be a positive radon measure on M, and denote by 



L„(X) = J exp(Ax)/i(eb) (1.1) 

its Laplace transform. Let Ai(JR) be the set of measures [i such that the set 

6(/i) = interior{A G JR.; L M (A) < +00} (1.2) 

is non empty and [i is not Dirac measure. The cumulant function of an element \x of 
A4(JR) is the function defined for A 6 0(/x) by 

U\) = hxLJ\). 
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To each \i in A4(IR) and A in 6(/x), we associate the probability distribution on IR 

P(A, /i)(dx) = exp(Ax - fc M (A))/i(dx). (1.3) 

The set 

F(/i) = {P(A,/i); AeeW} 
is called the natural exponential family (NEF) generated by \x. 

The function /c^ is strictly convex and analytic. Its first derivative k' defines a diffeomor- 
phism between 0(/x) and its image M F ^y Since fc^(A) = J x u )(dx), M F ^ is called 
the domain of the means of F(fi). The inverse function of k'^ is denoted by tp/j, and setting 

P(m,F( M )) = P(^(m), M ) (1-4) 
the probability of with mean m, we have 

= {P(m,F(,t))i mG M F(jU) }, 

which is the parametrization of by the mean. 

The variance of P(m,F(/j,)) is denoted by Vp( M )(m) and the map defined from M F ^ into 
L S (H), the set of symmetric function in R, by 

m i — > V F(p) (m) = k'^ip^m)) = (^(m)) -1 

is called the variance function of the NEF F(/j,) generated by /x. We also say that a 
is a basis of F(fj,). An important feature of Vpv^) is that it characterizes the natural 
exponential family F(fi) in the following sense: If F(fi) and F(y) are two NEF's such 
that V F ^{m) and V F ^(m) coincide on a nonempty open subset of M F ^ n Mp^j, then 
F{v) = F{v). 

Now, for |U 6 A4(R) the J0rgensen set of /Lt or of is defined by 

A(jjl) = {A > 0; 3 L^(9) = (L^8)) x and 6(/x A ) = 9( M )}. 

A(/i) is stable under addition which means that if A and A' are in A(/z), then A + A' are in 
A(/x), and u x+y = a x * /U A '- 
For all A in A(/x) we have 

TTl 

M F(nx) = XM F(,t) and V F(llx) (m) = W FM (—). 

Several classifications of NEFs based on the form of the variance function have been 
realized in the last three decades. The most interesting classes of real NEF's are the 
class of quadratic NEFs, i.e., the class of NEF's such that the variance function is a 
polynomial in the mean of degree less than or equal to 2 characterized by Morris [11], 
and the class of cubic NEF's, i.e., the class of NEF's such that the variance function 
is a polynomial in the mean of degree less than or equal to 3 characterized by Letac 
and Mora [10]. Recall that up to affine transformations and power of convolution the 
class of quadratic NEF's contains six families: the gaussian, the Poisson, the gamma, 
the binomial, the negative binomial, an the hyperbolic family. The class of cubic NEF's 
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contains, besides the quadratic ones, six other families, the most famous is the inverse 
Gaussian distribution with variance function V(m) = to 3 . It is worth mentioning here 
that multivariate versions of these classes have also been defined and completely described. 
For instance, Casalis [1] has described the so-called class of multivariate simple quadratic 
NEFs and Hassairi [6] has described the class of multivariate simple cubic NEFs which are 
respectively the generalizations of the real quadratic and real cubic NEF's. The fact that 
the variance function of a family is quadratic or cubic is not only a question of form but 
it corresponds to some interesting analytical characteristic properties. Indeed, the Morris 
class of quadratic NEF's has some characterizations involving orthogonal polynomials due 
to Fiensilver [5] . These characterizations have been extended to the Letac and Mora class 
of real cubic NEF's by Hassairi and Zarai [8] using a notion of 2-orthogonality of a sequence 
of polynomials. Other remarkable characterizations of the quadratic NEF's are related to 
the Bayesian theory. For instance, given a NEF F(p), Diaconis and Ylvisaker [4] have 
considered the standard family II of priors on the natural parameter A defined by 

(dX) = C tl , mi exp(tiTOiA - hk^X)) 1©(^)(A) dX (1.5) 

where t\ > 0, toi is in M F ^, and Ct 1>mi is a normalizing constant. This distribution is in 
fact a particular case of the so called implicit distribution on the parameter of a statistical 
model introduced in [7]. They have shown that if AT is a random variable distributed 
according to P(X, fi) (see (1.3)), then the only conjugate family of prior distributions on A 
that gives a linear posterior expectation of k'^(X) given X is the standard one II (see also 
[2]). Consonni and Veronese [3] have considered another family II of prior distributions 
7ft! , mi on the mean parameter m defined also for t\ > and toi in M F ^ by 

TT tl , mi (dm) = C tumi exp(tiTOi^(TO) - hk^^m))) 1 M {m)dm. (1.6) 

F{u.) 

They have shown that the fact that II contains k'^H) characterizes the quadratic NEFs. 
These authors have also shown that if the prior on the mean parameter to, is 7r tl , mi , then 
under some conditions on the support of /U, the NEF F{(i) is quadratic if and only if the 
posterior expectation of k'^X) is a linear function of the sample mean. We also mention 
that Diaconis and Yilvisaker [4] have shown that if the standard prior on A is given by 
TTti.mi with t\ > and vri\ is in M F ^, then the expectation of A;^(A) is equal to mi, that 
is 

C tl ,mi J k'^(X) exp(£imiA — hk^X)) l & ^(X)dX = mi, (1.7) 
or equivalently in terms of the mean parameter 

/TO, 
— -— exp^imiV^m) - tik^Jrn))) 1 M (m)dm=m 1 . (1.8) 

A natural question within this approach is if one can extend the properties and character- 
ization results concerning the quadratic NEF's and related to the Bayesian theory to the 
Letac-Mora class of real cubic NEFs. The aim of the present paper is to give an answer 
to this question. We first introduce, for a given NEF F(v) and (3 in some interval of IR 
containing 0, a NEF F^(u) such that F°(u) = F(u). We then define a family H 13 of prior 
distributions on the natural parameter 9 and a family YlP of prior distributions on the 
mean parameter to which may be seen as generalizations of the families II and n defined 
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above, since II = II and II = II . After proving that for each /3, the family is a 
conjugate family of prior distributions with respect to the NEF F"(u), we show that a 
cubic NEF F{v) is characterized by the fact that there exists a j3 such that the posterior 

expectation of tf, \„ s is linear when the prior on 6 is 7rf „ . We also show that a cubic 

1 - (3k' v {0) ' 

NEF F{y) is characterized by a differential equation verified by the cumulant function 
k v . A third characterization of a real cubic NEF is realized when the family of priors 
contains the family k' u ijl^). The restriction of all these results to the subclass of quadratic 
NEF's leads to the results of Diconis and Ylvisaker[4] and Consonni and Veronese [3]. The 
results of the paper are illustrated by an example. 



2 Main results 

In this section, we state and prove our main results. Our considerations will be restricted 
to regular NEFs, so that the domain of the means of a NEF is equal to the interior of 
the convex hull of its support. This property of regularity is satisfied by all the most 
common NEF's. An important fact which will be crucial in our proofs is that up to affine 
transformations and powers of convolution, a cubic natural exponential family may be ob- 
tained from a quadratic one by the so-called action of the linear group GL(\R 2 ) on the real 
families. Originally, the action of the linear group includes the affine transformations and 
powers of convolution, however since these transformations preserve the class of quadratic 
NEF's and the class of cubic NEF's, we will focus on the facts which we need here, for 
more precise statements in this connection, refer to Hassairi[6]. Let F(y) and F(p,) be two 
real NEFS's. Suppose that there exists a /3 in H such that the set 

{M F{v) )p = {m G Mp{y) ; 1 + /3m > 0} 

is nonempty and for m in (Mp)p, 

V F(v) {m) = (1 + /3m) 3 V m {jT^j ' (2 ' 9) 

then we write F{y) = Tp(F([i)). This defines an action on the natural exponential families, 
so that we have TpTpi = Tp + pi and F(y) = Tp(F((i)) is equivalent to F(fx) = T_p(F(u)). 
We also mention that F{y) = Tp(F([i)) may be expressed in terms of the cumulant 
functions of the generating measures by 

' K{9) = fc„(A) 

(2.10) 

6 = -/3^(A) + A 

or equivalently by 



(2.11) 



A 



,(0) + 



An important fact is that the relation (2.11) between the cumulent functions may be 
explicitly given it terms of the measures themselves. In fact if the a-power of convolution 
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v a of v is written as v a (dx) = h(a, x)a(dx), where a(dx) is either Lebesgue measure or a 
counting measure, then the measure 

»(dx) = - _ 1 ^ x /t ( 1 - /^^AH^ - (3x)a(dx) 

satisfies (2.10) and generates the family T-p(F(v)). This measure fj, will be denoted T-p{v) 
and the family = T_p(F(v)) will be denoted F' 3 . 

We mention here that if F(y) is a cubic NEF, there exists (3 in .Bp^) and a quadratic NEF 

such that F(u) = F(Tp{n)) or equivalently, F(u) = F(T_p(u)). 
Besides the restriction to half lines for the domain of the means in the definition of 
(M F t v \)p, we also define for v G M(M) and ft G M, the sets 

= {x G IR; 1 + /3x G A(i/)}, 

and 

B F{y) = {/3 G M; > 0}. 

We have the following preliminary result. 

Proposition 2.1 Let F(v) be a regular NEF. Then 

ft G B F{u ) if and only if (M F ^)p / 0. 

Proof We will make a reasoning for /3 > 0, the case /3 < may be done in a similar 
way. Suppose that there exists mo in (M F ^)p, that is too G M F i v \ and 1 + /3m-o > 0. As 
M F i v \ is equal to the interior of the convex hull of supp(u), there exist xo in supp(v) such 
that xo > m,Q. This with the fact that 1 + /3too > imply that 1 + /3xo > 0. Thus Hp is 
an open set which contains an element of supp(ji). It follows that v{Hp) > and (3 is in 

B F (y). 

Conversely, if /3 is in B F ^, then v{Hp) > 0. Since iTg is an open set, this implies that 
it contains an element xo of supp{v). We have on the one hand that 1 + (3xq > so that 
there exists e > such that 1 + /3xo — fie > 0. On the other hand, as M F ^ is equal to the 
interior of the convex hull of supp(u), there exists mo in (Mp(„j) such that |mo — xo| < e. 
From this we deduce that mo is in (M F ^)p. □ 

We now use the natural parametrization and the parametrization by the mean of the origi- 
nal family F(y) to give two parameterizations of the family F(ji). These parameterizations 
are, for j3 / 0, different of the usual parameterizations of F(fi). In fact, for 8 G 0(V), we 
write 

P(/3, 0, v){dx) = exp{(0 + /3£v(#))x - fc„(0)} T^p{u){dx). 
Similarly, parameterizing by m G M F ^ , we write 

P(0,m,F(i/))(dx) =exp{(^(m) + /3^(^(m)))x-^(^(m))} T_0(i/)(dx). 

Thus we have that 

F(ji) = F? = {P(/3, 6, u)(dx); 9 G 0(i/)} = {P(/3, m, F(!/))(dx); m G M F(l/) }. 
Accordingly, we define for /3 in B F t u \ two families of prior distributions. Let 

{Q)p = {9e@(u); l + [3kUe)>0}. 
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Then we have that (Mp^)^ = ^((0)^), and we define for t > and tuq G (Mp^p, 

4„(«) = Cfmo (1 + ^(5)) exp(im o - ifc„(0)) 1(0),^)^, (2.12) 

and 

^ = {^i t > and m G (M F(iy) ) /3 }, 

This family comes in fact from the standard family II defined in (1.5) using (2.11). The 
normalizing constant Cf mo is then well defined for t > and mo G {M F r v \)p. 
Besides this family of priors on the parameter 9, we define a family of priors on the pa- 
rameter m. Always for t > and mo G (Mp^)^, we consider the probability distribution 

i„W = ^5n„ (l + /3my 2 exp(tmo^„(m) - tk^^m))) l (AfpM)/j (ro)dro, 

where C^ m() is a normalizing constant. It is the image of 7Tti,mi defined in (1.6) by the 
m' 

map m' i — > — -. The family of priors on m is then 

1 — pm! 

& = {^t, mo \ i > and mo G (M F(i/) ) /3 }, 

Next, we prove that these families are conjugate families of prior distributions. 

Proposition 2.2 i) The family IP is conjugate family of prior distributions with respect 
to the NEF F 13 parameterized by the natural parameter 9. 

ii) The family IP is a conjugate family of prior distributions with respect to the NEF 
F 13 parameterized by the mean parameter m. 

Proof i) Suppose that X is a random variable with distribution P((3,9, v)(dx) and that 
7r^ mo is a prior distribution on the parameter 9. Then the posterior distribution is 

Cf mo (1 + PK(e)) exp((tm + x)9 - (t + 1 - gaQA^g)) l (&h (9) 
J Cfmo (1 + Pk'M) exp((im + x)9 - (t + 1 - f3x)k^{9)) l (e)j9 (0) d9 

tmo ~\~ x 

If we set ti = t + 1 — (3x and m,2 = — , then this distribution is noting but 

t + 1 — px 

^t+i-px (tmo+x)/(t+i-px)> an d ^ belongs to II 13 . In fact, since t > 0, T^p{v) is concen- 
trated on {1 — (3x > 0}, and mo is in (M F ^)p, we have that ti > and 1 + /3m,2 = 

1 + t(l + /3m ) , N 

— — — — — > 0, so that m 2 is m {Mf)^. 

t + 1 — pA 

ii) Suppose now that Y is a random variable P((3,m, F(v)) distributed and that the prior 
on the mean parameter m is 7f^ mo , then the posterior distribution of m is given by 

Cl mo (1 + /3m)" 2 exp((im + y)^(m) - (t + 1 - py)k^{m))) \ Mv{v) )^ m ) 
j Ct,mo (i + fimy 2 exp((tmo + i/)^(m)- (t + 1 - py)k^^m))) l ( M F(l/) ) /3 ( m ) dm 

This with t 3 = t + 1 -/3y and m 3 = -^jL±Jj- is equa i to Srf^^^+yy^^, with 
the required conditions to belong to IP\ □ 
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Corollary 2.3 Let (X±, X n ) be a sample P (ft, 9, is) -distributed and consider ^t,m as 
the prior on 9. Then the posterior distribution of 9 given X\, X n is 

t+n-f3nX,(tm +nX)/(t+n-l3nX) 

Proof It is easy to see that the distribution of the random vector (9, Xi, X n ) is 

n 

Cfm (1 + PKW) exp((tm + nx)9 - (t + n(l - px))K(9)) l (0)/3 (0) J[ T_p(v)(dxi)d9. 

i=l 

With the same technic used in Proposition 2.2, we deduce that the posterior distribution 
of 9 given X\, ...,X n is 



7I "* +n- /3nX , (tm +nX ) / (t +n-/3nX ) ' 



□ 



Proposition 2.4 Let -F(^) 6e a cubic NEF. Then there exists f3 in B F ^ such that, if the 

L - ; K(P) \ _ "'0 



prior on 9 is vr^ mo , then 



l + f3k' u {6)J 1 + /W 



Proof Since F(v) is cubic then there exist j3 in B F t v \ and a quadratic NEF F(fi) such 
that F(z^) = Tp(F(n)) and f = Tp(n). Using (2.11) it is easy to see that if the prior on 
in (0)/3 is ^t,m then the prior of A in 0(/x) is the standard 7i"t limi with ii = i(l + /?mo) 
and mi = - ^ ■ Moreover we have Cf m() = C t(1+/ 3 mo))mo/(1+/3mo) . 
It follow that 

^(irik) = /^r^w (1+ ^ WW(w -^ ))1(e) ' Wd9 

= / <n ^(A) exp(tm (A - ^(A)) - %(A)) 1 0(m) (A)(2A 
= y C tumi k^(X) exp(timiA - ti^(A)) l e(ju) (A)dA. 



Invoking (1.7) we get 



E 



k' v {9) \ m 



l + Pk' v {e)J 1 + /W 

□ 



Now we give a characterization of the cubic NEFs which is based on the linearity of the 
posterior expectation. 

Theorem 2.5 Let v be in M(M). 

1. If F(y) is cubic then there exists j3 in B F f u \ such that for all n > 1, if X±, X n 
is a sample with distribution P(/3, 9, v) and the prior on the natural parameter 9 is 
Kt, mo , ^en E I Xx, AT n ) is linear. 
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2. The converse is true if we assume that supp{T-p{v)) contains an open interval in 
M or is a finite or denumerably infinite subset of] — oo,0[ ; alternatively ]0, +00 [, 
with T-p(i/){0} > 0. 

Proof 

1. From Proposition 2.4 and Corollary 2.3, we easily get 

E ( Kid) ^ \ = tm + nX 



1 + PK(0) ' n ) t{l+fim )+n 



which is linear in X. 



2. Conversely, suppose that there exists fi in B F ^ such that E ( \^pk\e) I ^l, ■■■■>X n 
is linear for all n and a sample X±, ...,X n with distribution P(fi, 6, u)(dx). Consider 
fi = T-p(y). Then using (2.11), we have that the distribution P(fi,9,u)(dx) of 
Xi,...,X n is equal to P(X, fi)(dx) which is nothing but an other parametrization 
involving fj,. We also have that 

which is linear in X±, ...,X n , from the hypothesis. 

On the other hand, as the prior on the natural parameter 6 is assumed to be 7r^ mo , 
we get as prior on A G the standard prior distribution given by 

n u mi(d\) = C tl , mi exp(timiA - tife M (A)) l & ^(X)(dX), 

with ti = t{\ + fim ) > 0, mi = - ^o mQ G M F{^) and c h, mi = Cf mo - 

As we have that supp(fi) = supp(T^p(v)) C supp(v) n {1 — fix G A(i/)}. the as- 
sumptions on supp{T-p(v)) imply that supp(/j.) satisfies hypotheses (HI) or (H2) 
of Theorem 1.1 of Consonni and Veronese. According to this and to the linearity 
of the conditional expectation of the mean parameter of F(fi), we deduce that this 
NEF is a quadratic. It follows that F{v) = T (3 (F(^)) is a cubic NEF. 

□ 

In the following theorem, we give a second characterization of the Letac-Mora class of real 
cubic NEFs. 

Theorem 2.6 Let v be in M.(JR). F(y) is cubic if and only if there exist fi in B F ^ and 
(a,b,c) G -K 3 such that for all m in Mp^j, 

V F{u) {m) = (1 + fimf exp (aip u (m) + &fc„(^„(m)) + c). (2.13) 

Note that (2.13) may be expressed in terms of the cumulant function as there exist fi and 
(a, b, c) G R 3 such that for all in 0^ 

k'l(0) = (1 + fik'Mf exp (a8 + bk u (9) + c), 
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that is the cumulant function is solution of some Monge- Ampere equation (see[12]). Proof 
Suppose that F{y) is cubic, then there exist /3 in B F r v \ and a quadratic NEF F(ji) such 
that F{v) = Tp(F(n)) or equivalently Fiji) = T_p(F(v)). Then it follows from (2.9) that 



V F{y) {m) = {l+PmYV F ^[- 



m 



(3m 



It is known (see [1]) that for the quadratic NEF F(n) there exists (a',b',c r ) £ IR 3 such 
that for all m! in M F ^ 

Writing (2.11) in terms of the mean parameters we get 

' K (W(ttM = k v {Mm)) 

< 



(2.14) 



Therefore 

Vf{v) ( m ) = (1 + /5m) 3 exp(aip u (m) + bk u (ip v (m)) + c), 

with a = a', b = b' + /3a' and c = c'. 
Conversely, if (2.13) holds, then 

In V F ( u )(m) = 31n(l + /3m) + aip v {m) + bk u (^ u (m)) + c. 

Taking the derivative, we deduce that the variance function satisfies the differential equa- 
tion 

(1 + l3m)V F{v) {m) - 3pV F(u) (m) = (a + bm){l + /3m). 
Solving this equation by standard methods gives 

V F{y) {m) = A(l + /3m) 3 + A(l + /3m) 2 -^^(l + /3m), 

which is a polynomial of degree less than or equal to 3. □ 

A third characterization of the cubic NEF's is based on a relation between the associ- 
ated families of prior distributions IF 3 and IF. 

Theorem 2.7 Let v be in Ai(M). F(v) is cubic if and only if there exist j3 in B F ^ such 
that 

k' u {T{P) C &. 

Proof Suppose that F(u) is cubic, then from Theorem 2.6 there exist /3 in B F ^ and 
(a, b, c) in IR 3 such that 

V F{u) (m) = (1 + /3m) 3 exp (a^(m) + bk u (ip u (m)) + c). (2.15) 

Consider the set 

n = {(ti,mi) € Hi; x {M F{ff) ) p -M - b > and € (M F(i/) )^}. 
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We will show that k' v {Ii p ) = {vff ljmi ; {h , mi) G Jl}, which is a part of IT 3 . 

Let t > and mo in (M^Vj,))^, and denote by a the image by k' u of the prior 7r 3 m(J on 
defined in (2.12). We easily verify that 

a(dm) = Cf mo (1 + /3m) V~^ u) {m) exp(tm Vv(m) - tk v (il) v {m))) l^^^mjdm. 

This using (2.15) becomes 

cr(dm) = Cf i mi (1 + /3m)' 2 exp{timiip u (m) - hk^tp^m))) l^^^imjdm, 

with U = t + b, mi = — - — - — and Cf ™, = Cf , /. ^ ,„\/u m- 

1 t + 6 { i,mi ti~ b\timi+a) / (t-i — b) 

We have that t\ — b = t > and - = mo G (M f i v \)r, that is (ti,mi) G fi. Hence 

t\ — b 

^(n^)c{7ff limi ; (ti,mi)eO). 

In the same, we verify that 

(ti,mi)€fi}c*an"). 

Finally, we obtain that 

= i^mi, (h,m 1 )en}c& 

Conversely, suppose that k' v {nP) C II 13 . The image of an element 7rf mo of IT' 3 by fcj, is by 
the very definition 

Ki.^t,m )( dm ) = c t, mo (1+>S?tc) (Vj^m)) -1 exp(im Vv(rn)-iMVv(m))) 1 (M FM ) (9 ( m ) dm - 
Since it is assumed to be in IF, there exists (t\, mi) in JR* + x (M F i v \)p such that 

K{^im ){dm) = C? umi (l + /3m)" 2 exp(iimi^(m) - hk^ip^m))) l {Mv)fj (m)dm. 
Comparing these two expressions of k' v (vr^ mo ) gives 

V F(u)( m ) = (1 + /3m) 3 exp(aVv(m) + bk u (ip u (m)) + c), 

where 

/ r-/3 \ 

/ t Tna i 

a = tmo — tiirii, b = t\ — t, and c = In' ' 1 



According to Theorem 2.6, this is the desired result and the proof is complete. □ 
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3 Example 



In this section we illustrate our results by an example involving the most famous family 
with variance function of degree 3 which is the inverse Gaussian natural exponential family. 
Consider the distribution 

(1 + x) _3//2 x 2 

v(dx) = == exp( - -) lii , nn r(x)dx : 

v J ^ v\ 2 {l + x)' 1 1,+oo[V ' ' 

which is up to an affine transformation an inverse Gaussian distribution. The NEF gen- 
erated by v is given by 

F{y) = {exp((l + x)9 + V^29) u(dx); 9 < 0}, (3.16) 

its mean parametrization is 

n») = + rhi ] m > (3 - 17) 

For all m > —1 the variance function is given by 

V F{u) (m) = (l+m) 3 . (3.18) 

Let \i = T-p{y), we have 

" {dx) = V^^W+- X eXp( -2(l + ( X l-/3)x) ) W-flW*)** 
Now for f3 7^ we have 

p0-\-y/—26 ™2 

P{ ^ 6 ^ dx) = V^Vl-Px + x eM{e+(3{ ~ e ~ VZ ^ ))X ~ 2(l-^x + x) ) Hl-^O}^, 
so that 

= T^(F(u)) = {P(0, 6, u)(dx); 9 < 0}, 
The corresponding family II 13 of conjugate prior distributions is the family of distributions 

4 mo (d8) = Cfmo (1 + ^7=ff " *)) ^P{tm 9 + t(9 + ^29)} l (&h (9)d0, 

where 

(©)/3 = ]-^(^l) 2 '°[' 



]-l//3,+oo[ i/ /3<0 

]-l,+oo[ if (3 = 

]m/(-l,-l//?),+oo[ i/ /3>0 



and 

(CU) = "J [ X - ex P(-^-(^3T) +^3T-2(^3T) )]-(l+^mo) |_^ ( _^_ )2 - 
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defined for t > and mo in (Mpf v \)p 

Also, in this example, the family IF is the set of distributions defined for t\ > and mi 

in (M F(i/) )^ by 

To see how Theorem 2.6 holds in this example, we need only to take (3 = 1. Then jjl = 
T_i(i/) is the standard gaussian distribution with Vp^(m!) = 1 for m! G H. We see that 

^F(i/)("i) = (1 + m ) 3 exp(a^(m) + bk v (ifj v (m)) + c), 

with a = b = c = 0. In fact 

v F(ji){™') = exp(aV M (m') + b'k^tp^m')) + d), 

with a' = b' = d = 0, and using the relations a = a', b = b' + (3a', c = d, we get 
a = b = c = 0. 

Concerning Theorem 2.7, we fist observe that the hypotheses in this theorem are well 
verified. In fact, let 7r^ mo be the prior on the natural parameter 9 and 7^ mi be the prior 
on the mean parameter m, then 

tj = t + 6 = t > 
tmo — a 

mi = ^ = m £ (Mf (i/) )i =J1,+oo[. 

The density function of ^t,m * s e Q ua l to 

* exp{im 6> + v/^)}, 

v — 20 

and for all m > — 1 the density function of k' v {ir\ mo ) is given by 

/-, \-2 r / I 1m 

(1 + m) expj ; -77 — t(— -77 H 

V ; P1 2(1 +m) 2 V 2(l+m) 2 1 + m" 

which is equal to nl mo . 

References 

[1] Casalis, M. (1996). The 2d+4 simple natrual exponential families on IR d , 
Ann. Statist. 24 1828-1854. 

[2] Cifarelli, D.M. and Regazzini, E. (1983). Qualche Osservazione sull'Uso di dis- 
tribuzioni iniziali coniugate alia famiglia esponenziale, Statistica, 43, 415-424. Cor- 
rections: Statistica(l990), 50, 293. 

[3] Consonni, G. and Veronese, P. (1992). Conjugate priors for exponential families hav- 
ing quadratic variance functions, 

Journal of the American Statistical Association, 87, 1123-1127. 



12 



[4] Diconis, P and Ylvisaker, D. (1979). Conjugate priors for exponential families. 
Ann. Statist 7 269-281. 

[5] Feinsilver, P.(1986). Some classes of orthogonal polynomials associated with martin- 
gales, Proc. Amer. Math. Soc. 98 298-302. 

[6] Hassairi, A. (1992). La classification des families exponentielles naturelles sur ]R rf par 
Taction du groupe linaire de H n+1 , C.R.Acad.Sc. 315. 

[7] Hassairi, A. Kokonendji, C, and Masmoudi, A. (2004). Implicit distribution and 
estimation, Communication in Statistics, 34 2, 245-252. 

[8] Hassairi, A. and Zarai, M. (2004). Characterization of the cubic exponential families 
by orthogonality of polynomials, Annals of probability, 32 2463-2476. 

[9] Letac, G. Lectures on natural exponential families and their variance function. Mono- 
grafias de Matematica. 50, IMPA, Rio de Janeiro. 

[10] Letac, G. and Mora, M. (1990). Natural real exponential familiies with cubic variance 
functions, Ann. Statist, 18 1-37. 

[11] Morris, C. N. (1982). Natrual exponential families with qudratic variance functions. 
Ann. Statist. 10 65-80. 

[12] Zuily, C. (1989). Local existence and regularity of the Dirichlet problem for the 
Monge- Ampere equation, Journes "equations aux Drives Partielles" (Saint Jean de 
Monts, 1989), Exp. No. XXI, 5 pp. 



13 



